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Abstract 

I summarize Density Functional Theory (DFT) in a language familiar to quantum 
field theorists, and introduce several apparently novel ideas for constructing systematic 
approximations for the density functional. I also note that, at least within the large I\ 
approximation (K is the number of electron spin components), it is easier to compute 
the quantum effective action of the Coulomb photon field, which is related to the 
density functional by algebraic manipulations in momentum space. 
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1 Density Functional Theory for Quantum Field The¬ 
orists 

Much of atomic, molecular and (quantum) condensed matter physics, reduces, in the non- 
relativistic limit, to the problem of solving the Schrodinger equation for point-like nuclei and 
electrons, interacting via the Coulomb potential. The electrons are fermions, and are best 
treated by introducing an electron held 

^t(x) = ^a;(k)e ik ' x , ( L1 ) 

k 


where i is a spin label. If we introduce dimensionless space-time coordinates, by using the 
Bohr radius to measure length, and the Rydberg to measure energy, the only parameters in 
the problem are the nuclear charges Z a and the mass ratios The latter appear multiplying 
the nuclear kinetic terms 


\ - m e 

11 kin—nuc / . 17 

^ m a 


( 1 . 2 ) 


a 


The mass ratios are all < .0005, and this leads to the Born-Oppenheimer approximation 
for calculating a systematic series in powers of these small dimensionless parameters. One 
first freezes the nuclear positions and finds the electronic ground state energy VBo(R-a) as a 
function of the nuclear positions. Then one minimizes the Born Oppenheimer potential and 
expands around the minimum to quadratic order. 

The minimum value of Vbo is the order 1 contribution to the ground state energy of the 
system in Rydbergs. The oscillator corrections are down by s /— . 

The next stage depends on whether one is dealing with atoms, molecules or large crys¬ 
tals. For atoms, the ground state is spherically symmetric. For molecules, the minimum 
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of the Born Oppenheimer potential has a shape (the one you find models of in chemistry 
departments), and can be rotationally excited. The rotational level splittings are of order 
—. For large crystals the coherent superpositions of rotational levels corresponding to fixed 
positions of the crystal, satisfy almost classical equations of motion. On the other hand, 
the low energy spectrum of vibrational levels reaches down to \fm e /m a /L, where L is the 
linear size of the crystal. These phonon excitations are important to a number of properties 
of solids. 

The hard quantum mechanical part of all of these problems is the solution of the electronic 
Schrodinger equation for fixed nuclear positions. It is a many fermion problem with no 
dimensionless parameters. For infinite crystals one can introduce the electron density in 
Bohr units as a dimensionless parameter. Density Functional Theory is a general approach 
to this difficult problem. It was invented by Kolin Hohenberg and Shams [lj. 

In this section, I want to present an elementary approach to DFT using concepts familiar 
to particle theorists, in the hope that we can make some contributions that may have been 
missed in the condensed matter literature. DFT has had remarkable successes in all areas 
of atomic, molecular and condensed matter physics, but much of the work is numerical. My 
hope is that there are some more analytical approaches to the calculation of the density 
functional, which can be of some use. In the second part of this paper I will present an 
analytic approach, based on an expansion in the number of electron spin components. 

The Born Oppenheimer approximation reduces all of our problems to that of electrons 
interacting via Coulomb repulsion, in a background potential V ext , which neutralizes their 
total charge. V ex t is a finite sum of nuclear Coulomb potentials, but it is convenient to 
generalize it to be an arbitrary function of x . Our problem then is to find the ground state 
energy of interacting electrons in a general neutralizing external potential, as a functional of 
the potential. 

Field theorists will instantly recognize this as the problem of finding the generating 
functional kF[14//(x)] of connected equal time Green’s functions of the electron density 
operator 

W(x) = V’iWV’iW, (1-3) 

in the Homogeneous Electron Gas (HEG). This observation is, from one point of view, 
the foundation of density functional theory. The density functional is just the functional 
Legendre transform of W. 

More intuitively, we can use the variational principle to fold the ground state energy by 
minimizing 

(x\ J d 3 x'0l( x )(-V 2 -/r)V’i(x) + ^ J d 3 x d 3 y + J d 3 x N(x)V ext (x),\x), (1-4) 

over all normalized states |x). is the chemical potential and 7V(x) = 7V(x) — n 0 , where no 
is the neutralizing charge density. We perform the minimization in stages, first minimizing 
the expectation value of the Hamiltonian among normalized states with a fixed value of the 
density 

™( x ) = (x\N(x)\x), (1.5) 

and then minimizing the resulting functional of n(x) w.r.t. the density. 
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This procedure recommends itself, because the non-universal part of the energy, the part 
that depends on the external potential, is just 



and is thus a simple and explicit functional of the density alone. The hard part of the 
calculation, which gives rise to the density functional F[n(x)] is a property of the HEG: it is 
the minimum expectation value of the HEG Hamiltonian in states that have fixed expectation 
value of the density operator iV(x). 

The message is clear, if surprising. The better we can do at solving the apparently aca¬ 
demic problem of the HEG, the better we can do at all the problems of atomic, molecular, and 
condensed matter physics. Since the ground-breaking Hohenberg-Kohn paper that pointed 
this out [2] condensed matter physicists and quantum chemists have been making this insight 
pay off in a major way. The HEG itself is a difficult problem, with no intrinsic dimensionless 
expansion parameter besides the average electron density in Bohr units. For high density, an 
expansion around the free electron gas is useful, though it must be re-summed to deal with 
apparent IR divergences. The point is that the long distance Coulomb potential is screened 
(the quantum analog of Debye screening in a classical charged gas at finite temperature) so 
the long distance and weak coupling limits do not commute. This problem is easily dealt 
with and several terms in the high density expansion of the ultra-local part of F[n] are 
known. 

Much of the work on DFT uses an approximation Erst proposed by Kohn and Shams 
(KS) [3]. They write 

F[n) = F 0 [n} + ^ f d 3 x d 3 y M x Wy) + p LDA ^ (1.7) 

^ J l x y| 

where the first term is the density functional of a free electron ga^j], the second, called 
the Hartree term, is the interaction term of the Hamiltonian with the density operator 
replaced by its expectation value, and the third, which is usually calculated by sophisticated 
numerical algorithms ra. is an ultra-local functional of the density, ft is the ground state 
energy density of the HEG. A variety of fitting functions have been proposed, which fit the 
numerical data and match the high and low density expansions HU, but there does not seem 
to be a systematic approach to computing corrections. There have been various attempts 
to go beyond the local density approximation. Some particularly successful ones are found 
in p2], but the most accurate fit to real data is found by making empirical fits to parameters 
in theoretically motivated formulae [T3J. There is, to my knowledge, no systematic expansion 
for the density functional. One aim of the present paper is to propose directions for such 
a systematic set of approximations. We will see that the leading term in an expansion in 
the inverse number, K , of electron spin components gives something very close to the first 
two terms in the KS ansatz. The only difference is that F 0 [n] is expressed in terms of the 
density, rather than KS orbitals. However, at the density which gives minimum energy, the 

Where is a subtlety in this statement, in that Kohn and Shams write Fq [n\ in terms of single particle 
solutions to the Schrodinger equation — V 2 + /z— = Ci'bo rather than the density itself. We’ll deal with 

this subtlety in the section on the KS equations. 
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two approaches give the same result. The higher order terms in the 1/K expansion are all 
describable in terms of diagrams whose vertices are functional derivatives of the partition 
function of a non-interacting electron gas w.r.t. a time dependent background potential. The 
propagator is similarly related to the density two point function. These remarks suggest an 
important step forward will be to develop analytic functional approximations to the non¬ 
interacting problems with an external held. 1 will make some tentative remarks about this 
in the penultimate section of this paper. 

The large K expansion may also make it possible to explore analytically the properties 
of low density states of the HEG, which violate translation symmetry spontaneously, like 
the Wigner crystal [3], and emulsions or stripes [5j. The leading order large K expansion 
gives a steepest descent approximation to the expectation value of the potential coupled to 
the density operator. This is a translation invariant equation, but can have solutions which 
violate translation invariance and minimize the energy. Large K expansions are not always 
reliable guides to the detailed physics of such phases and the transitions between them, but 
they certainly provide a convenient starting point for quantitative discussion. We’ll see that 
in leading order, the large K expansion predicts a first order transition between a Wigner 
crystal and a homogeneous phase, at zero temperature. 

The KS approximation certainly gets some qualitative features of molecules and solids 
right. If we combine the Hartree term, the external potential term, and the Coulomb repul¬ 
sion between nuclei, and ignore everything else, then the energy is minimized by making the 
electron density exactly cancel the nuclear charge density. However, because of the uncer¬ 
tainty principle, Fq is large if the density varies on scales much shorter than the Bohr radius, 
while the nuclear wave functions are localized on scales shorter by a factor of \Jm e /m a , so 
exact cancelation is impossible. However, the Bohr radius of a single electron ion with a 
nucleus of charge Z is 1/Z in Bohr units, so the electron density tends to cluster in the 
vicinity of nuclear positions and to try to screen out the large nuclear charges, over distances 
of order 1 in Bohr radii. Once the nuclear charges are screened down to something of order 
1, the competition between attraction to the nuclei, mutual repulsion, and Fermi statistics, 
becomes less easy to estimate, and one must go to more refined approximations. However, 
it’s clear that the qualitative fact that only a few valence electrons per nucleus are responsible 
for molecular binding, and transport in solids, follows simply from this crude approximation. 

What’s especially remarkable is that this essentially classical minimization problem cap¬ 
tures so much of the quantum mechanics of a many electron system. Field theorists are used 
to the platitude that, in principle, everything one wants to calculate can be “reduced” to the 
problem of finding stationary points of the quantum effective action, but in atomic, molec¬ 
ular, and condensed matter physics, these fantastical claims actually lead to manageable 
computations. 

To conclude this introduction, 1 want to address myself to quantum chemists, band 
structure theorists, or other professional utilizers of DFT, who might happen across it. This 
paper is not addressed to you. As a practical tool, the methods proposed in this paper are 
not ready for prime time. My main motivation in writing has been to get professional held 
theorists interested in the held theory of the HEG, which is enormously simpler than the 
relativistic theories we typically study, but of enormous practical utility. The professionals 
have developed DFT into a tool of impressive accuracy, but the foundations of the theory 
and a more systematic approach to it have not been a center of interest for a while, as far 
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as I am aware. A. Zangwill has pointed out a very interesting early paper [B], which takes 
a held theoretic approach to DFT, similar to that advocated here. It will be interesting 
to compare the approach taken there with my own. The application of modern techniques 
of QFT to the HEG, is certainly a great mathematical physics project, and it might lead 
to real advances in practical calculations. I’m taking my first tentative steps in this held, 
and I would certainly appreciate guidance to the existing literature and things I might have 
missed. 

Some textbooks/monographs on DFT can be found in [lTj, and an excellent recent history 
of the subject, from which I’ve benehtted enormously is [9]. 


2 Expansions in the Number of Spin Components 
2.1 The Large K Expansion 

From the point of view of Quantum Electrodynamics, the Hamiltonian of the HEG re¬ 
sults from integrating out the electro-magnetic held in the non-relativistic approximation for 
electron motion. We can re-introduce the scalar potential by a Hubbard-Stratonovich trans¬ 
formation. The finite temperature free energy of the HEG is then written as a Euclidean 
functional integral, with action for the Coulomb photon held 0 : 

S = f dr ^2 k 2 0k0-k + 2Tr In ( d T - V 2 - igcj)(T,x.) - V ext (x)). (2.1) 

■' k 

The trace is taken over the space of fermionic wave functions, -0(r, k) , periodic on a spatial 
torus and anti-periodic in Euclidean time. We’ve introduced the torus in order to isolate 
the spatial zero mode of the Coulomb photon held qb. This mode is omitted from the sum, 
corresponding to the fact the the zero mode of the nucleon charge density, exactly cancels 
that of the electron charge density for a neutral system. We absorb the chemical potential 
into V ex t . The dimensionless parameter g is equal to 1 when we work in Bohr-Rydberg units. 
We insert it in order to obtain a smooth large K limit when we give the fermions K spin 
components. 

The log of the fermion determinant, like most functional determinants, has an additive 
infinity. We normalize it by insisting that for r independent 0 it just equals the free energy 
of electrons in the external potential V — i<p. 

We now modify the HEG by allowing the electrons to have K spin components, turning 
the SU(2) spin rotation symmetry into SU(K). The density operator N(x ) = ip 1 ^(x)^i(x) 
has an expectation value of order K and connected two point correlations of order K , so, 
in order to have a smooth large K limit, we set g = yj2/K, and define 0 = a/ K/2o. The 
resulting functional integral over a has an action 


SH = it'd ^2 k^k^-k + Tr 111 (d T -V 2 -i(tj-iV at )] = /V j / a(-V 2 )cr + PL[cr-iV at ] 


( 2 , 2 ) 


2 This is done all over the condensed matter literature. I’ve followed some concise lecture notes by Ben 
Simons [7], with small changes in notation. 
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For large K the fluctuations of o around the stationary point of the action are of size 
-j=. For an equilibrium state, the stationary point should be independent of Euclidean time. 
Shifting variables to y = icr + V ext , the equation for a stationary point becomes 

V 2 (xo~V ext ) = --^-. (2.3) 

<5y(x) 

The functional derivative of L on the RHS is the electron density in a thermal equilib¬ 
rium state of non-interacting electrons with external potential yo- Note that in terms of 
the original variable a, the stationary point is on the imaginary axis, corresponding to a 
real potential. This often happens when evaluating an integral by the method of steepest 
descents. 

The value of the free energy E (defined by Z — e~ /3E ) at the stationary point is 


E[V ext ] = Ii 


d 3 x d 3 y 


- 5L 1 5L 

-<5y(x) —2V 2 X ' y 5y(y) 


+ L[xo\ , 


(2.4) 


where [3 is the inverse temperature. The functionals are evaluated at the stationary point, 
at real yo, which is independent of Euclidean time. 

The density n(x) is defined as the functional derivative of E w.r.t. U. To evaluate it we 
use the equation for the stationary point to write 


SU(x) 

$xo(y) 


= <5 3 (x — y) + / G(y,w) 


5 2 L 


^Xo(w)5y 0 (x)‘ 


We then calculate 


5E 

W) 


= / d 3 y 


5E 8xo( y) 
^Xo(y) SU(x) ' 


using the chain rule, and obtain 


6E 5L 

— K - 


SU(x) dxo(x) 
The free energy can now be written 


(2.5) 


( 2 . 6 ) 


(2.7) 


E{V „,] = / d 3 x cl 3 y | k (x ) _ 2A . V 2 (x. y)w(y) + KL [xo] 


( 2 . 8 ) 


where 


n(x) = K 


5L 

^Xo(x) ‘ 


Since KL[x o] is just the free energy of the system of non-interacting electrons in the external 
potential y 0 we can write 


Now use 


KL[x 0 ] = F 0 [n] + 


V 2 (y 0 - V ext ) = 


J [ n Xo] ■ 

(2.9) 

SL 

5y(x)' 

(2.10) 
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We see that the Legendre transform of E[V ext ], which defines the leading large K contribution 
to the density functional F[n] is just given by 


F[n] = F 0 [n] + ^ J d 3 x d 3 y ra(x)^(x, y)ra(y), (2.11) 

which is the standard Hartree approximation to the density functional if K = 2. Note that 
since n = 0 (K), both of these terms are of order K in the large K expansion. 

The term of order 1 in the large K expansion of the free energy is 


F x = -Tr In 
2 


-V 2 5 3 ( 


y)S(t x - ty) + 


6 2 L 


8x(x,t x )6x(y,t v 


( 2 . 12 ) 


The operator whose determinant is evaluated here is called the inverse propagator of the y 
field. 

Higher order terms in the large K expansion involve Feynman diagrams whose vertices 
come from higher functional derivatives of L, and propagators for the x field, which are 
inverses of the above integro-differential operator. The key to evaluating these corrections 
is the solution of the single electron problem in an external potential, as a functional of the 
potential. 

ft’s also interesting to understand the relationship of the density functional to the func¬ 
tional T[y], which field theorists call the effective action or 1PI (one particle irreducible) 
generating functional of the field y. This is defined in the HEG with no external potential. 
The diagrams of the 1/K expansion have a propagator that is given, to zeroth order by 
[—V 2 + n(x — a;')]” 1 , where I’ve introduced space-time notation to save typing (x° = t) and 
n is the one loop, Euclidean time dependent two point function of the density operator. 
One can divide all connected Feynman diagrams into those which cannot be cut in two, by 
cutting a single propagator line (the one photon irreducible diagrams), and those which can. 
It’s well known to field theorists that one can easily (meaning using only algebra in momen¬ 
tum space, without doing integrals) calculate the exact answer for any n - point correlation 
of the field y, if one knows the exact 1 PI correlators of y for k < n. For the two point 
function, the result is known as the Dyson series. One simply corrects, order by order in the 
1/K expansion, the polarization function Yl{x — x '), defined as the sum of all 1PI two point 
functions, and writes the full propagator as (—V 2 + n) _1 . 

1PI n-point functions, r n (xi ... x n ) of the field y are written in terms of Feynman graphs 
involving loops of one component free electrons, which define the k point vertices of the 
graphs as the connected, time dependent, k point correlators of a gas of free electrons with 
a single spin component. The k point vertex has a weight , The propagator for 

the photon contains only a single fermion loop contribution (RPA), and we sum over all 1PI 
graphs. The generating functional of 1PI graphs is 


rfc] = 


^ 1 


n =0 


n! 


d 4 nx r n (xi... x n )xc{xi) ■ ■ ■ Xc{x n ). 


(2.13) 


It is well known to quantum held theorists that 

r [Xc]+ J d 4 xJ(x)xc(x), 


(2.14) 
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evaluated at the value of \ c satisfying 


= (2.15) 

0Xc{x) 

is minus the logarithm of the functional integral over y with action 

= K \\ I xV 2 X + L[x} + J x{x)J{x). (2.16) 

The equation above says that T[y c ] is the Legendre transform of the negative log of this 
functional integral. Call the Legendre transform of T, W[J]. 

On the other hand, the density functional (generalized to time dependent density fields) 
is the Legendre transform, w.r.t. V ex t of the negative log of the functional integral with 
action 

S\x] = K{\ J(X- V ext )V 2 ( X ~ V ext ) + L[x\. (2.17) 

Thus F[n] is the Legendre transform of 

E[V] = W[-V 2 V ext ] + ^J ( V ext )V 2 (V ext ). (2.18) 

The Legendre transform is not a linear operation. However, the Legendre transformed n 
point functions, are given by tree diagrams constructed from the n point functions of the 
original functional (as vertices) connected by the inverse two point function. The functionals 
E and W have n > 3 point functions which are the same np to multiplication by powers of 
momenta of the external legs, and two point functions which differ by a shift by A; -2 .It’s 
easy to construct the momentum space n point functions of F[n] as algebraic combinations 
of the 1PI n point functions of T, thus exploiting the diagrammatic simplification of one 
photon irreducibility. 

An even easier procedure is simply to evaluate T[y c ] and use the relation 

E[V] = W[-VV e *t] + ^J ( V ext )V 2 (V ext ), (2.19) 

and the Legendre transform relation between W and T to evaluate E[V] directly, as the 
solution to a variational problem involving T. F[n] has the virtue of being manifestly gauge 
invariant under space independent, but time dependent, gauge transformations of y, while 
T is not. However, the simplicity of this calculational scheme suggests that it is the right 
way to proceed. 

2.2 Small K Expansion 

We can also contemplate a small K expansion of the functional integral over a. For this 
purpose it is convenient to go back to the original variable (f), in terms of which we see that the 
potential term in the Schrodinger operator becomes very large in the small K limit, while the 
whole determinant contribution becomes small. We know that neglecting the determinant 




entirely will not work for small wave numbers, where screening effects are important. This 
suggests that a good approximation to the determinant for large, slowly varying potentials 
will be useful in this limit. This is precisely where the WKB approximation to the fermion 
determinant should work. 

If indeed we succeed in constructing a small K expansion as well as a large K expansion, 
then the technique of “two point interpolation” (a generalization of two point Pade approxi- 
mants to expansions that are not pure integer power laws) might be expected to give reliable 
results for all values of K. 


2.3 The KS Equations 

All of the applications of DFT are based on a series of approximations to F[n] introduced by 
Kohn and Sham (KS) |3] . For non-interacting fermions, the density functional is computed 
by introducing independent single particle orbitals ^(x), in terms of which the density is 

n(x) = J]^*(x)V>, 0 (x). (2.20) 


These are the lowest N normalized eigenstates of the single particle Schrodinger equation 


- V 2 * + 11 +V (x)</>,° = £($. 


( 2 . 21 ) 


The value of N is taken to coincide with f n(x), which is tuned by tuning the chemical 
potential fi. The energy is 


J2 ei= [ X^*°)*( x )[~ V2 + h + E(x)]'0°(x) 

i i 


F 0 [n(x)] + 


V(x)n(x). 


( 2 . 22 ) 


From this we see that for the non-interacting problem 

FoN= / £>°)*(x) [-V 2 + ##?(*). (2.23) 

i 

Note that these -0*° are fixed functions, because the potential in the Schrodinger equation, 
which they satisfy, is determined by 

fly = ~ v(x) - {2M) 

Fq [n] is then simply the kinetic energy of these single particle orbitals, and the density is 
written n(x) = Ei(^°)*( x )0; ( x )- 

KS now take the single particle orbitals as their variational parameters, rather than the 
density. The density functional is taken to be 

F[n) = K[if* , ipi] + F Hartree [n\ + F xc [n ], (2.25) 

where n(x) = E«'0«*( x )f ; j( x )) an d K the kinetic energy functional of the orbitals. One then 
varies F[n] + f Vbo n w.r.t. the if* after making a suitable approximation to F xc , obtaining 
the KS equations: 

(—V + fjLWi - [ -- — - Vbo\A = (2.26) 


9 




The e* are Lagrange multipliers, which enforce the condition that the orbitals are normalized 
on the variational problem. 

These equations are solved by iteration. One starts with the solutions without the Hartree 
and X-C contributions to the density functional, computes the density, and plugs it into the 
density functional to get a new estimate of the self consistent potential, then repeats the 
exercise. My understanding is that the numerical convergence of the iteration is fast and 
stable. 

When satisfy the KS equations, the functional is not equal to the non¬ 

interacting density functional F () [n ], because the potential in the Schroodinger equation 
satisfied by the ijji is not equal to — ^. For general n there is no simple relation between 
the two functionals. However, when n satisfies the density functional condition 

SF 0 1 (F jj ar f ree [nj T F X C ) T T I 

M^) + [ -- Bo] ’ 


(2.27) 


then the KS orbitals indeed satisfy 

5F 

(—V 2 + n — ° )^i = e^i. (2.28) 

on(x) 

We conclude that, at the minimum, the KS functional coincides with that obtained by 
using F 0 [n] rather than the orbital kinetic energies. The KS equations are equivalent to 
those obtained by varying the density^. The KS procedure makes the term F 0 [n] a rather 
explicit functional of their variational orbitals, while computing F 0 [n] is a nontrivial project, 
because it’s tantamount to solving the single particle Schrodinger equation for an arbitrary 
potential. KS reserve this work for their numerical solution of the variational equations. 
Historically, this was surely motivated by the existence of large computer codes for solving 
the Hartree equations at the time KS began their work. 

In order to carry out computations in the large Ii expansion, we need to compute 
L[x(t, x)]. For time independent x this is the Legendre transform of F 0 . Thus, for our 
purposes, working in terms of n rather than KS orbitals, would seem to be the route to 
follow. 


3 Expansions of the Functional Determinant 

Functional determinants can be computed simply in two different limits, constant fields of 
arbitrary strength, leading to the derivative expansion, and weak fields of general functional 
form. We begin by showing that these two expansions are expansions in a single parameter 
around 0 and infinite values. Consider the Euclidean Schrodinger equation 

[d t - V 2 - n - gxit , x)]/ n = A n fn- (3.1) 

For small g we can solve this by simple perturbation theory. To understand the large g limit 
we rescale t —>■ t/g , n —)■ gg, x —> 1/gx. In terms of the new variables, the equation is 

[d t - v 2 - n - x(tg~\*g~ 1/2 )]fn = (3.2) 

3 I suspect that this is well known to adepts in DFT, but I could not find a clear explanation of this point 
in the literature. 
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The overall rescaling of the eigenvalues leads to an additive, x independent term in the 
logarithm of the determinant, which we can discard. It’s clear from this formula that the 
large g limit is equivalent to one in which the fields are slowly varying functions so the 
derivative expansion, implemented via the heat kernel expansion as in [6], is the appropriate 
approximation. The leading term gives the Thomas Fermi approximation to the density 
functional. 

ft’s certainly possible to compute several terms in the small g expansion fairly easily. The 
large g expansion is more complicated but we can certainly get the first non-leading term 
with relative ease. My goal is to find an interpolation between these two approximations 
which will be valid for all g , with small error. There are two different ways to think about 
such an interpolation. One can calculate physical quantities like the ground state energy 
of a particular system as expansions in g and 1/g and find interpolating formulae for each 
quantity individually. Alternatively, one can try to find an interpolating formula for the 
entire functional L[x], and then plug that into the 1 /K expansion. The second procedure is 
more universal, and adheres more closely to the original spirit of DFT, but it may be difficult 
to implement. 

4 1/K expansion of the HEG 

Let us now apply the large spin expansion to the ground state of the HEG. The leading 
order leads to a translation invariant equation for the background Coulomb photon field a 

2 1 SL 

° ~ T <5cf(x) ’ 

where the determinant is calculated with V ext = /i and a background a independent of 
Euclidean time. It is therefore proportional to the Euclidean time interval, T, which is taken 
to infinity. The tilde over a reminds us that the zero momentum mode is decoupled from the 
determinant, because of the constraint of neutrality. The only parameter in the equation is 
the chemical potential, or equivalently the expectation value of the density operator. 

The homogeneous solution of this equation, a = 0 gives the correct ground state for high 
density. On the other hand, an intuitive argument shows that for low density, a periodic 
solution, in which there are electron density fluctuations of order N above and below the 
neutralizing homogeneous positive background charge, is energetically favored. The scaling 
argument, which enabled us to scale out the Bohr Rydberg units, shows that the Hartree term 
in the equation is more important than the non-interacting electron density functional at low 
density. Intuitively, the Hartree term gives negative energy from a regular array of positive 
and negative charge excesses. One pays a price for this in electron kinetic energy because 
the band structure in the self consistent periodic potential gives a higher Fermi energy than 
entirely free electrons. However, as the length of the lattice spacing goes to infinity, this 
effect goes away, and the scaling argument shows that it goes away more rapidly than the 
negative Coulomb energy. 

Thus, at sufficiently low density there is a state with negative energy of order K, which 
is the true ground state, since the order K contribution to the homogeneous ground state 
energy vanishes. The quantum phase transition between these two systems is of first order. 


(4.1) 
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The inhomogeneous phase sets in at a finite wavelength, determined by the critical density. 
Gaussian fluctuations around the homogeneous phase are controlled by the Random Phase 
approximation to the Coulomb photon polarization. The operator 

p 2 + n (p,u), 

where II is give by a single fermion loop , has no negative modes for any value of the density, 
so the homogeneous phase is meta-stable in the large K limit. This is consistent with the 
finite, o(K), jump in energy between the two phases. 

I have not studied the details of the equation closely enough to do more than quote 
conventional expectations (see the Wikipedia article) for the precise form of the lattice that 
minimizes the energy at low density. These are the body centered cubic lattice in three 
dimensions, the triangular lattice in two, and a uniform lattice in one dimension. The large 
K limit gives a controlled context in which the arguments that lead to these expectations, 
should become rigorous. 

Another point, which one should be able to investigate quite simply in this limit, is 
whether there is a tricritical point in the density temperature plane, at which the transition 
becomes second order and disappears. Certainly, we expect that at high enough temperature, 
only a homogeneous phase exists. In the large K approximation, the critical quantity to 
study is the temperature dependent Coulomb photon inverse propagator p 2 + II (p,u>,T) in 
the homogeneous phase. As one approaches the critical point along any trajectory in the 
density-temperature plane, this operator should develop a zero mode and then an unstable 
mode as one passes through in directions where the crystal phase exists. It’s entirely possible 
that this computation already exists in the literature, so I will put it off to a future paper 
until I’ve done a more thorough literature search. 

Higher order terms in the large K expansion around the homogeneous phase are given in 
terms of Feynman diagrams. The propagator in these diagrams is the inverse of p 2 + n(p, u>), 
and the vertices are connected correlation functions of the density operator in the free electron 
gas. n itself is the two point connected density correlator in the free gas. In a future paper, 
I intend to compute a few orders in the large K expansion of the ground state energy of the 
HEG, and compare them to the accurate numerical answers from quantum Monte Carlo. 

5 Conclusions 

I’ve summarized DFT in a language familiar to quantum field theorists and proposed several 
strategies for finding analytic approximations for the density functional. I’ve also suggested 
that the effective action for the Coulomb photon field might be an even simpler quantity 
to calculate. It seems to me that there is a lot more that can be done, and that field 
theorists should be able to do a really accurate job of calculating at least the connected n 
point functions of the density operator of the HEG for n < 4 and the first few terms in the 
derivative expansion of the effective action, to all orders in the density. 

Ignorant of the fundamental results of DFT, I’d always considered the HEG an academic 
problem of limited interest, but I now realize how central it is to all ab initio calculations of 
the properties of matter under the normal circumstances that prevail on earth. It’s worth 
putting a lot of effort into this problem. I hope others will agree. 
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